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In this paper we analyze the delay in a single-server queue in 
which the server, when it becomes free, selects for the next service the 
oldest customer with current delay smaller than T. If no such cus- 
tomer is present, then it selects the youngest customer with the current 
delay in excess of T. This service discipline is desirable in applica- 
tions where the success or failure of a service depends on the delay in 
providing the service. Telephone call processing and steel rolling are 
two of these applications. We obtain the delay distribution for this 
service discipline using a combination of level-crossing arguments 
and renewal theory, and compare this performance with that of the 
last-in-first-out discipline with respect to the throughput of success- 
fully served customers. 

I. INTRODUCTION 

The following situation is common in telephone call processing or 
data-processing systems. When a customer requests service, an entry 
is made in the queue that is serviced by a processor. The processor 
serves the queue of entries according to some specified service disci- 
pline. When an entry is served, the corresponding customer is notified 
of the completion of the service. The customer, however, does not wait 
forever for the completion of the service. At some random time, R, 
after its arrival, the customer will renege if service is not completed. 
The associated entry remains in the queue, and the server does not 
know that the customer has reneged until after it completes the service 
and attempts to notify the customer. Such a service is wasted. The 
customer may make the situation worse by reattempts at getting the 
desired service, thereby increasing the load on the server. It is therefore 
necessary to keep the proportion of reneging customers as small as 
possible. This can be done by selecting an appropriate discipline. 

In this paper we study a specific model of the above situation. In 
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Fig. 1 — Customer behavior function 1. 

particular, we have a single-server queue with Poisson arrivals at rate 
A>0. 

Let G be the distribution of the service time. Let P(t) be the 
probability that the customer does not renege before time t. We may 
think of P(t) as the expected reward obtained by completing the 
service of a customer at t time units after its arrival. With this general 
interpretation P(t) need not be restricted to be between and 1. Let 



Pit) 



-I 



P(t + y)dG(y) 



(1) 



for < t < oo. Then P(t) is the expected reward from a customer whose 
waiting time (excluding the service time) is t. Let W„ be the distribu- 
tion of the waiting time under the specified service discipline, m. Then 
the expected reward from an arbitrary customer is 






P(t)dWAt). 



(2) 



We want to select a service discipline that maximizes V„. It was shown 
in Doshi and Lipper 1 that if P(t) is convex (respectively, concave), 
then the last-in-first-out (LIFO) [respectively, first-in-first-out (FIFO)] 
discipline is optimal. More realistic functions P(t), however, are of the 
forms given in Figs. 1 and 2. 

For such functions, P(t), an optimal service discipline, is not known. 
However, our results for concave and convex P(t) indicate that a 
hybrid discipline may provide better performance than either the 
FIFO or the LIFO discipline does. In this hybrid discipline the server, 
when it completes a service, first looks at the customers with the 
current waiting time less than Tand selects the oldest waiting customer 
for the next service. If no such customer is waiting, then the server 
looks at the customers with the current waiting time in excess of T 
and selects the youngest customer. Note that this hybrid discipline 
includes FIFO (T = oo) and LIFO (T = 0) as special cases. Since P(t) 
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Fig. 2 — Customer behavior function 2. 

is assumed to be given, we only need W„(t), < t < oo, for this hybrid 
discipline it, to calculate V„. In this paper we obtain W„ for general 
service time distribution. We do this in three steps. First we describe 
another queueing system for which the distribution of the waiting time 
is the same as for the original system. Moreover, for this equivalent 
queueing system the distribution of the waiting time can easily be 
expressed in terms of the distribution, F(x), of the work in a subsystem. 
We then use level-crossing arguments to derive an integral equation 
satisfied by f(x) = F'(x). Finally, we use some results from renewal 
theory to solve this integral equation. 

Some comments about the model are in order before we proceed to 
give an outline of the rest of this paper. Models similar to the one 
studied here can be useful in a variety of other applications. Some of 
these are the management of steel-rolling operation and the manage- 
ment of blood bank. Also, in many applications the customers do not 
necessarily renege. They simply take actions (start to dial, become 
cold, etc.) which make any subsequent service worthless. 

This paper is organized as follows: In Section II we formally define 
the queueing system under consideration. We describe an equivalent 
queueing system in Section III. There we also show the relationship 
between the distributions of the waiting time in the original system 
and of the work in a subsystem of the equivalent queueing system. In 
Section IV we derive an integral equation for the steady-state density 
of the work in the subsystem. We give the solution of this integral 
equation in Section V. There we also derive the steady-state distribu- 
tion of the waiting time in the original system. Finally, we give some 
numerical results in Section VI. 

II. MODEL 

The queueing system and the hybrid service discipline discussed in 
Section I can be formally described as follows: We have a queueing 
system with a single server and two queues, Ql and Q2. Customers in 
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Ql have a nonpreemptive priority over those in Q2. The order of 
service is first-in-first-out (FIFO) in Ql and last-in-first-out (LIFO) in 
Q2. Customers arrive according to a Poisson process at rate X > 0. On 
arrival the customer is put in Ql. If its service has not started within 
T seconds after its arrival, then at that time it is transferred to Q2. 
The service times of the customers are independent and identically 
distributed with distribution function G with continuous density g. 

Let fix denote the Kth moment of the service time and let p = Xfii . 
Assume that the waiting-time process is in the steady state. Let W 
denote the distribution function of the waiting time seen by an arbi- 
trary customer. Since G has a continuous density, W is differentiable 
on (0, oo). Let 

w(x) = W'(x) < x < oo. 

We are interested in obtaining an expression for W{x) or, equivalently, 
for W(0) and w(x), < x < oo. 

III. AN EQUIVALENT QUEUEING SYSTEM 

We now describe a queueing system that is equivalent to the one 
described in Section II as far as the waiting times of the customers are 
concerned. However, the number of customers in Ql and Q2 at a given 
time may be different in the two systems. 

Consider the subsystem consisting of the server and Ql. Let X t 
denote the work, at time t, in this subsystem. Thus, X t is the sum of 
the remaining service time of the customer, if any, being served and 
the service times of all the customers in Ql. If a customer arriving at 
time t finds X t ^ T, then it joins Ql; otherwise it joins Q2. Recall that 
Ql has a nonpreemptive priority over Q2 and that Ql is served FIFO 
and Q2 is served LIFO. A little reflection shows that the waiting time 
of a customer is the same in this system as in the one described in 
Section II. 

We now relate the waiting-time distribution, W, to the steady-state 
distribution, F, of the work, X, in the subsystem consisting of the 
server and Ql. If an arriving customer sees X < T, then its waiting 
time will be X because the service in Ql is FIFO and because Ql has 
a nonpreemptive priority over Q2. Thus, 

W(x) - F(x) < x < T. (3) 

In particular, 

W(0) = F(0), (4) 

and 

w{x) = f(x) = F'(x) 0<x<T. (5) 
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If an arriving customer sees X > T, then it joins Q2. Since Ql has 
priority over Q2 and since Q2 is served LIFO, this customer has to 
wait for the duration of a busy period in the M/G/l queue with initial 
work, X. Let B y denote the distribution of the busy period in the 
M/G/l queue started by initial work, y. Since G has a continuous 
density, 



exists for all x > y. Then 



by{ X ) = B'y(x) 



w(x)= f(y)b y (x)d(y) x>T. 
Thus, it is sufficient to find the distribution of X 



(6) 



(7) 



IV. INTEGRAL EQUATION FOR f (x) 

We now use level-crossing arguments (see Ref. 2) to derive an 
integral equation for f. Figure 3 shows a typical sample function for 
the process {X-}. Assume that at time the queues are empty, the 
server is idle, and a customer arrives and enters service. If X t > 0, then 
it decreases at unit rate as in a M/G/l queue. Symbols o represent 
arrivals which see X t < T and join the subsystem, thus increasing X t 
by a service time. A symbol in the shape of a dot (•) represents 
arrivals that see X, > T, and join Q2 without affecting X, on their 
arrival. When X t reaches zero, two things can happen: Q2 is empty and 
the server remains idle until the next arrival, or Q2 is nonempty and 
a customer from Q2 enters service, thus increasing X t by its service 
time. Such arrivals from Q2 into the server are denoted by a symbol in 
the shape of a square (■). 

The stochastic process {X t } is not Markovian because what happens 




Fig. 3 — Typical sample function for the process {X}. 
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when X, reaches zero depends on the past. On the other hand, the 
vector- valued process {(X t ,N t )}, where Nt is the number in Q2 at time 
t, is a Markov process. We can derive the steady-state distribution of 
(X t , N t ) using the standard results and use that to obtain the marginal 
steady-state distribution of X. We, however, use a simpler approach 
here. Recall that /X2 < °°- Consider the following two cases: 

(i) p < 1. In this case the process {X t } is regenerative and with the 
regeneration points corresponding to the external arrivals that make 
an idle server busy. Denote this event by E. Then E is a positive 
recurrent, regenerative event. 

(ii) p > 1. In this case the queue length in Q2 grows without bound 
and, in the steady state, can be assumed to be oo. Thus, a customer is 
removed from Q2 to enter service every time X t reaches 0. This event, 
E', is then positive recurrent. 

Standard regenerative arguments now show that {X} has a steady- 
state distribution and that 

Xt-^X, 

where the distribution of X is the steady-state distribution of {Xt) . 
Moreover, for any x, < x < oo, the steady-state rate, D(x), at which 
Xt crosses x from above, equals the rate U(x), at which X t jumps from 
below x to above x. We now express D(x) and U(x) in terms of f(x) 
and get the desired integral equation by equating these expressions. 

X t decreases at unit rate until an arrival occurs or until X t = 0. Thus, 
during every downcrossing of level x, the {Xt) process spends dx units 
of time in the interval (x — dx, x). Hence, 

D(x) = f(x) < x < oo. (8) 

Before deriving an expression for U(x) we introduce some notation. 
Letp denote the rate at which Xt jumps from to some positive value. 
These jumps may be due to either external arrivals coming to an idle 
system or to customers from Q2 moving to the server. Also, let G 
denote the complementary service time distribution defined by 

G(x) = 1 - G(x) < x < oo. (9) 

Assume that G(0) = 0, G(0) = 1. Then, for x < T 



U(x) = A f(y)G(x - y)dy + pG(x). (10) 

Jo 

Since an external arrival causes a jump in X t only when X t < T, we 
have, for x > T, 

U(x) = A J f(y)0(x - y)dy + pG(x). (11) 

1256 THE BELL SYSTEM TECHNICAL JOURNAL, MAY-JUNE 1983 



Thus, 



Let 



•TAx 

U(x)=\\ f(y)G(x - y)dy + pG(x) 0<*<oo. (12) 



Jo 



f(0)=hmf(x). (13) 

x[0 

Then f(Q) is the rate at which {X t } hits level zero. In the steady state, 
this must equal the rate at which {X} jumps from to some positive 
value. Thus 

P=f(0). (14) 

We now have for 0<p<oo J 0<.x<oo, 

fix) = D(x) 
= U(x) 



Let 



= A f(y)G(x-y)dy + f(0)G(x). (15) 



/(0) 

9 = 1T- 



Then 



•xAT 

fix) = X I f(y)G(x - y)dy + q\G(x). (16) 



(*) - A J 
Jo 



This is the desired integral equation for f. The additional conditions 
needed to solve this completely depend on whether p < 1 or p > 1. 

First consider the case where p < 1. Let P2 be the probability that 
an arriving customer sees X, > T and joins Q2. Since, for p < 1, every 
arriving customer is eventually served, the rate at which customers 
enter the server from Q2 is XP2 . Also, the rate of arrivals coming to an 
empty system is XF(0). Thus, 

g Jf.M^±M. Fm + Pl . (1 7) 



Also, 



■I 



P 2 = f(x)dx, (18) 
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and 

F(0) + f(x)dx = 1. (19) 

Jo 

For p > 1, Q2 is always nonempty in steady state. Hence F(0) = 0, 

,-*?. <20) 

and 

/(x)rfx = 1. (21) 



Equation (16), together with either conditions [eqs. (17) to (19) or 
20], characterizes f completely. We solve this equation in the next 
section. 

V. SOLUTION OF THE INTEGRAL EQUATION 

We now solve eq. (15) to obtain an expression for fix). For < x < 
oo, let 

h(x) = XG(x) (22) 

and let m(x) be the renewal density function for h(x). Then mix) 
satisfies (see Ref. 3): 

m(x) = h(x) + I h(x- y)m(y)dy. (23) 

Jo 

Equation (15) can now be rewritten as 

fix) = qhix) + J fix - y)hiy)dy (0 < x < T), (24) 

Jo 

and 

fix) = qhix) + [ fiy)hix - y)dy T < x < oo. (25) 

Jo 

Equation (24) is a renewal equation and its solution is given by 



fix) = qhix) + q \ hix - y)miy)dy 



= Q 



hix) + hix - y)miy)dy 

Jo 

= qmix), (26) 
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where the last equality follows from (23). Note that < h{x) < w, and 

h(x)dx = Xfii = p. 



Thus, h is like a probability density with total mass, p. The function 
m(x) is well defined for any finite x irrespective of the value of p, < 
p < oo. To obtain an expression for f(x), x > T, we note that the right- 
hand side of (25) involves f(y) for y only in the interval (0, T), which 
we have obtained in (26). Thus, replacing f{y) on the right-hand side 
of (25) by qm(y), we get 

•T 



f(x) = 



= Q 



| qm(y)h(x-y)dy\ 
I h(x - y)m(y)dy 



h(x) + \ h(x- y)m(y)dy x>T. (27) 



We now use conditions (17) through (19) or (20) through (21) to 
evaluate q and thus completely characterize /. First, consider the case 
p < 1. We have 

q = F(0) + P 2 , (28) 

and 

• T 

F(0) + q I m(x)dx + P 2 = 1. (29) 



Also, equating the rate of customers coming to the system with the 
rate of customers leaving the system, we get 

X=-[l-F(0)l 
Mi 

or 

F(0) = 1 - p. (30) 

We can now solve (28) through (30) for q and P 2 to get 

f 
p — (1 — p) m(x)dx 

Jo 

Po = 



J-T 
m{x)dx 
o 



QUEUEING 1259 



and 



q = F(0)+P 2 -zt 



1 + 



r 

m(x)dx 
Jo 



Thus, 



F(0) = 1 - p 

,. . m(x) 

f(x) = -pf — 



1 + I m(x)dx 



I 



x<T 



f 

Jy=0 



h(x)+ h(x - y)m(y)dy 



1 + I m(x)dx 

Jo 

Next, consider the case where p>l. Here, 

/(0) 



9 = 



A ' 



and 



Jo Jt 



f(x)dx + f(x)dx = 1. 



Thus, 



(31) 
(32) 



(33) 



/■(0) | m(x)dx + p 1 + I m(*)dx - m(x)rfx | = 1, 



or 



/(0) = 



1 + m(x)dx 



f(x) = 



m(x) 



p 1 + 



f 

I m(x)rfx 

Jo 



*<T, 



(34) 



(35) 
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and 



/(*)- 



h(x) + h(x - y)m(y)dy 



1 + m(x)dx 



x>T. 



(36) 



We now consider a special case where G(x) = e l ",n= 1/jtii. Then, 



h(x) = Ae""* < x < oo, 

m(x) = Xe*-^* < x < oo, 

1+ I ro(x)dx = — !-— [1 - pe-" (1 - p)7 l 



I 



(1-P) 



and, for x > T, 



h(x)+ h(x - y)m(y)dy 



-Xe'^ + X 2 e^'-^e-^-^dy 



Jo 
= Ae-'" +Ar . 



Thus, for p < 1 



F(0) = 1 - p 



/(*)« 



1-P 



1 _ De -^T 



,-»(l-p)x 



x<T, 



and 



/(*)- 



1-p 



1 _ pe -M<i-p)^ 



,-jtx+AT 



*>T. 



For p > 1 



F(0) = 0, 
fix) = p[l - pe" ,p - 1)7 '] ' 



/iip-n* 



x< T, 
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and 



f(x) = 



1-P 
Ptl-pe"'"-^] 



,-fur+XT 



x> T. 



VI. WAITING-TIME DISTRIBUTION 

We now use the results of Section V to obtain an expression for the 
waiting-time distribution and its Laplace Stieltjes Transform. From 
eqs. (4), (5) and (6) we get 

1-p p<l 
.0 p>l, 

m(x) 



W(0) = F(0) = 



(37) 



w(x) = f(x) = -l 




P<1 



P2=l, 



(38) 



and 



h(y) + I h(y — z)m{z)dz 



b y (x)dy 



w(x) =< 



1 + I m(x)dx 
'o 

h(y) + I h(y — z)m(z)dz 
'o 



p<l 



(39) 



b y (x)dy 



1 + J m(x)<fc 

'0 



P>1. 



Let 



f 



W*(0) = E[e~ ew ] = W(0) + j e- 6x w{x)dx Re >0. (40) 
Also, for Re B > let 



£ 



Jo 



: g{x)dx 



(41) 
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h*(6) = e- 0x h(x)dx 



e 



(42) 



and 



m *{0)= e- Bx m(x)dx 



h*(0) 



l-h*{0) 

\[i-g*m 

6-\[l-g*(0)] 



(43) 



Let B*(9) be the Laplace Stieltjes Transform of the ordinary busy 
period in an M/G/l queue with the arrival rate X > and the service 
distribution, G. Then B * satisfies 

B*{0) =g*{6 + X[l - B*{6)]). (44) 

Also, let B$ (0) be the Laplace Stieltjes Transform of B y . Then 

B *( B) = e -v(fl+A[l-B-(fl)]) (45) 

We now express W*(6) in terms of m, m* and B*. First, consider the 
case p < 1, where 



e~ 8x m(x)dx 



W*(0) = l-p+ -pr 



1 + | m(x)dx 

10 



+ 



M 



= l-p + 



Jy=T L 

1 



h(y) + | M v - 2)m(z)d2 
e~ e *m(x)dx 



1 + m(x)rf^; 

'0 



9 
6 + \[l-B*(d)] 



b y (x)dyd. 



4 



e - (fl+ Ml-B-««»].-r m(;c)rf;c 



)■ 



(46) 
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For p > 1, 

•T 



W*(0) = — 4 { I e- 9x m(x)dx 



1 + I m(x)dx ' 



o _ r^.-B.^^j (47) 



+ \[l-B 

6.f Mean value of the waiting time 

Let W denote the mean value of the waiting time. For p < 1, every 
customer is eventually served. Hence, W is the average over all the 
customers. For p > 1, some of the arriving customers do not get served 
and, in this case, W is the average waiting time of those who do get 
served. In the first case, 

W=-W*'(0 + ), 

and for the second case, 

-W*'(0 + ) 



W = 



W*(0 + ) " 



For p < 1, all customers are served and the order of service does not 
affect the mean waiting time. Thus, in this case the mean waiting time 
is the same as that for an M/G/l queue with the FIFO discipline. That 

For p > 1, the busy period distribution is defective. Let 
b = P{Busy period < oo}. 
Then b D is the unique solution in (0, 1) of 

B*(0 + ) = b o = g*[Ml-b o )l (49) 



Also, 



and 



1 l + \g"[Ml - &„)] ( ' 



a»m = g '" [M1 ~ b ° )] (sit 

B (0> a + te-wi - wr ( ' 
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Now, from eq. (47), we get 

1 



W*(0 + ) = 



H 



m(x)dx 



m(x)dx + 



A(l - b ) e-o* 



lim Be- {0+X[i - B ' m])x m(x)dx 



r rT 



H 



m(x)dx 



m(x)dx + 1 





_1 
P 


and 




-W*'(0 + ) = 



(52) 



H 



m(x)dx 



xm(x)dx + 



Ml - &o) 



1 + m(x)dx 



}) 





" r 


1 


J xm(x)rfx 

Jo 


p 


f T 


1 + m(x)dx 




L Jo 



+ 



Ml - 6o) 



(53) 



From (52) and (53) we have 



W = 



xm(x)dx 



+ 



-W*W)__ _ 

W(0 + ) f r Ad - b„) 

1 + m(x)cfo 



(54) 



Equation (54) shows that the mean waiting time of the customers who 
get served is minimized by the LIFO discipline (T = 0). 

VII. NUMERICAL RESULTS 

In this section we present some numerical results. Instead of calcu- 
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lating the waiting-time distribution, we calculate the quantity of inter- 
est, namely 

V= I P(t)dW(t) 



-I 



P(t) = 



(55) 



for a specific P. Let 

t< T 

e ° {t - T) t > T. 

We now evaluate V for the hybrid discipline with parameter T and 
also for the LIFO discipline. For the hybrid discipline we have 

•T 

Vr=d-\ \ m(x)dx + e oT m*{o + \[l-B*(o)]} 



.(1 + h*{a + A[l - S*(ff)]» 

- e oT (l + h*{a + A[l - B*(a)]}) 

• J m(y)e- y(a+X)[1 - B ' {o)i dy 
Jo 



-U 



where 



1 + I m(x)dx 



Q = 



H. 



m(x)dx 



P <1 



P >1. 



(56) 



(57) 



(58) 



We need some more notation before writing an expression for Vl. 
Let 6fr( ■ ) denote the density function for the busy period started by 
the forward recurrence time of the service time. Thus, 



bvu(x) - I M*)- — — dy. 
Mi 



(59) 



Then, for p < 1, 



V L = 1 - P + P bv R (x) + pe° T e ° x b F R(x)dx 
'o Jt 



=l-p+p 



+ e ( 



r'V rT 

bFR(x)dx - e° T e~ ax b F R(x)dx 



A[l - B*(o)] 



o + \[l-B*(o)Y 
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(60) 



and, for p > 1, 

V L - I 6p R (*)rfx - e° T I c- x 6FR(*)rfx 

Jo Jo 

l-B*(o) 
Ml {a + A[l-fi*(a)]}' 



(61) 



For numerical calculations we considered all integrals in eqs. (56) 
through (58) and (60) through (61) as functions of T, obtained their 
Laplace transforms, and inverted the transform at the specified value 
of T using the inversion method of D. Jagerman. 3 Thus, let 

Ri(t) = m(x)dx, 

Jo 

R 2 (t) = I m(x)e- xlo+X[1 - B ' { °™dx, 



Rs(t) = brn(x)dx, 

Jo 

and 

R 4 (t) = I e-° x b FR (x)dx. 
Jo 

Also, for i = 1, 2, 3, 4, and in the appropriate domain, let 
R* (6) = I e- e 'Ri(t)dt. 



1 



Then 



„* (m m*{8+o + \[l-B*(o)]) 
K 2 (d) = 



and 



RU6) = 



d{6 + \[l-B*(6)]}' 

l-B*(6 + o) 

o + 6 + \[l-B*(9+o)] 



For numerical examples we had the service time distribution gamma 
with mean 1 and variance 1/K. We used two different values of K, 
K = 1 (exponential distribution), and K = 10. We used two different 
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values of T, 1 and 3. Finally, we used two values of a, 2.0 and 0.15. 
These give us eight parameter sets. The values of Vt and Vl as 
functions of the load A = p are given in Figs. 4 through 11. From these 
figures we observe: 

(i) For both service disciplines, the throughput of good calls is 
larger for larger T, smaller a, and larger K. The behavior of the 
throughput with respect to T and a is obvious. Larger K implies smaller 
variability in the service time, thus reducing the probability of a 
customer getting served after a long wait. This, in turn, results in a 
higher throughput. 

(ii) For the assumed customer behavior, the hybrid discipline al- 
ways provides higher throughput than the LIFO discipline does. The 
difference is larger for larger T, larger a, and larger K. 




Fig. 4 — The values Vt and Vl as functions of the load X = p for T = 3.0, a = 2.0, and 
K=l. 




Fig. 5 — The values Vt and V L as functions of the load A = p for T = 3.0, a = 2.0, and 
K = 10. 
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Fig. 6 — The values V T and V t . as functions of the load \ = p for T = 3.0, a = 0.15, and 
K=\. 




Fig. 7— The values V T and V; as functions of the load A = p for T = 3.0, o = 0.15, and 
K= 10. 




Fig 8— The values V T and V L as functions of the load \ = p for T = 1.0, a = 2.0, and 
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Fig. 9 — The values Vt and Vl as functions of the load A = p for T = 1.0, o = 2.0, and 
tf = 10. 




Fig. 10 — The values Vt and Vt as functions of the load A = p for T = 1.0, a = 0.15, and 
JP-1. 



Of course, our knowledge of the customer behavior may be more or 
less accurate, depending on the application. An issue of interest then 
is the sensitivity of the throughput to the assumed customer behavior. 
This was studied for a special case (K = 1) in Ref. 1. The analysis in 
this paper can be used to answer such issues for more general service 
time distributions. Qualitatively, however, the conclusions will remain 
the same: the last-in-first-out (LIFO) discipline is robust with respect 
to the knowledge of customer behavior. The hybrid discipline, on the 
other hand, is very sensitive to the customer behavior and should be 
used only when the customer behavior is adequately known and does 
not change in time, or when the parameters of the customer behavior 
can be estimated and used to change the control parameters dynami- 
cally. 
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Fig. 11 — The values V T and V L as functions of the load A = p for T = 1.0, a = 0.15, and 

K= 10. 
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